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Abstract 

We discuss the role and the treatment of polarization effects in many-body 
systems of charged conducting clusters and apply this to the statistical frag- 
mentation of Na-clusters. We see a first order microcanonical phase transition 
in the fragmentation of -Aa 70 + for Z = to 8. We can distinguish two frag- 
mentation phases, namely evaporation of large particles from a large residue 



and a complete decay into small fragments only. Charging the cluster shifts 
the transition to lower excitation energies and forces the transition to disap- 
pear for charges higher than Z = 8. At very high charges the fragmentation 
phase transition no longer occurs because the cluster Coulomb-explodes into 
small fragments even at excitation energy e* = 0. 

1 Introduction: statistical picture 

Atomic clusters are especially interesting as a link between macroscopic matter 
and a single atom. In studying clusters one would usually like to understand how 
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the material properties evolve with the size of the system. Our special interest 
is slightly different. We would like to address the very general question of the 
thermodynamics of finite systems. In this context finite means that the size of the 
system is of the same order of magnitude as the range of the forces. This imposes 
an overall correlation within the system which strongly influences its properties. 
Atomic clusters together with nuclei and gravitational systems are the examples of 
finite systems. 

This paper is the third in a series on "Fragmentation Phase Transitions in Atomic 
Clusters". In the first paper subtitled " Microcanonical Thermodynamics" i"l| we 
addressed the most general concepts, and in the second "Symmetry of Fission of 
Charged Metal Clusters" Q we work out how the material properties influence the 
fragmentation behavior. In the current work we investigate the influence of charge 
and polarization effects on the fragmentation of cold and hot Na-clusters and discuss 
it in light of microcanonical fragmentation phase transitions. The next work planned 
is on the relation of the fragmentation transition to the commonly known liquid-gas 
phase transition. 

Our overall goal is to calculate the accessible all-possible particle phase space Q 
after the fragmentation of a cluster within a given volume. The size of this accessible 
phase space is restricted by the conservation of the global parameters: mass, charge, 
total energy and the total angular momentum, which is set to zero in all calculations 
presented here. In addition there are geometrical constraints: the fragments within 
the cluster configuration should not overlap after the decay and, moreover, should 
not be closer than the binding distance of the daughter clusters. In studying a 
microcanonical system and applying Microcanonical Thermodynamics, we need to 
define the total volume which is accessible to the fragmenting system. The size of 
this fixed volume is also one of the geometrical constraints. 

In using the MMMC (Microcanonical Metropolis Monte Carlo) model []], H we 
simulate the accessible phase space at the moment after the fragmentation when the 
short-range Van-der-Waals forces (and the exchange of atoms) become negligible. 
This happens at the average surface distances between the fragments of 0.5 to 1 A. 



We simulate the fragmenting configurations within a spherical volume of radius 

1/3 ° 

Rsys — r f * Nrp , where N? is the number of atoms in the cluster and rf = 3.85A 
for Na-clusters. R sys is called the freeze-out radius. We want to mention in advance 
that the results of our calculations are practically independent of R sys , if r/ is kept 
within reasonable limits. We place all the fragments randomly inside this freeze- 
out volume under the constraint that the surfaces are further apart than a binding 
distance db ~ 0.5 to 1A (see Appendix C). 

It is also worth noting how the material dependence is included. The ground 
state binding energies of all possible fragments make up the input to our calcula- 
tion. For sodium we take the known experimental binding energies for neutral and 
singly charged clusters up to mass 21 |L8|, [7|]. Since we do not know the binding 
energies of other clusters we determine them from the liquid drop model. Thus the 
electronic shell effects are included only for smaller clusters, which are, however, the 
most relevant to the fragmentation. The weights of individual events (fragmenta- 
tion configurations) are also dependent on parameters like the Wigner-Seitz radius, 
moments of inertia of the fragments, and the vibrational frequencies of dimers and 
trimers. We treat the internal excitations for all the fragments with Nj > 4 as in the 
bulk matter where we use the bulk specific entropy (see references |j], |2|] for details). 

In studying the Coulomb explosion it is important to take polarization effects into 
account. We do this by calculating the pairwise interactions between the fragments, 
which we assume to be spherical (See Appendices A and B). The monomers are 
considered not to be polarizable. 

2 Coulomb explosion 

The investigation of the stability of charged water droplets goes back to the work of 
Lord Rayleigh |§ in the year 1882. His ideas were widely applied to similar prob- 
lems in nuclear and atomic cluster physics. Afterwards, scientists in both of these 
fields were interested especially in binary fission induced by Coulomb instability as 
a possible decay channel |§. In our calculations we allow all possible decay chan- 



nels of highly charged clusters. This enables us to study the competition between 
evaporation, fission, multifragmentation, or even total vaporization. 

The first experimental evidence for Coulomb instability of clusters was presented 
by Sattler et al., [|I| for small doubly charged clusters. Some experimental tech- 
niques for achieving highly charged states for masses N > 1000 are discussed in 
| II]. Here we want to concentrate on masses less than 100 in order to study highly 



charged states by including the polarization effects and their connection to general 
questions of thermodynamics in small systems. A comparison with experimental 
results is a matter for the future. 

Generally, the microcanonical phase transition in a small system is signaled by 
an S-shape in the caloric equation of state T(E) || [12], |13] , or in other words, a 



decrease in temperature with growing energy of the system in the transition region: 
the system cools while it absorbs energy. The thermodynamic temperature T is 
obtained by calculating the specific entropy, s, from a knowledge of the phase space 
Q(ttot), where e to t = E/N is the specific total energy: 

"4 - £ « 

ln(Q(etat)) 

s = ^v— 

For a fixed mass and charge of the cluster, f3 depends only on the specific excitation 
energy e*, 

^tot — C + £binding{N, Z) 
F N,Z 
( AT r Z\ — binding / \ 

CbindingK 1 * , ^ ) ~ TT • {*) 

e to t is the total specific energy available to the fragmenting system and is a globally 
conserved quantity, e* is also a conserved quantity , but only for fixed N and Z. In 
this case, considering the temperature dependence, it is sufficient to look at T(e*). 

We study the fragmentation of clusters at a fixed volume, as we assume that the 
motion of complex fragments at short relative distances is dissipative and ergodic. 

In order to concentrate on Coulomb and polarization effects, we have examined 
the statistical fragmentation behavior of the cluster Na 70 + at different charges Z. 



Figure |l| shows the fragmentation of Na^ for a calculation done assuming pairwise 
interactions between fragments represented by conducting spheres for a minimal 
surface distance d& = 1A (See appendices). The left scale shows the three largest 
masses Ml, M2, M3 of the fragments (averaged over all configurations) as a function 
of the specific excitation energy e*. We see that at e* below 0.2 eV the system mostly 
evaporates monomers. As e* increases, some heavier particles start to appear, but 
the most dramatic effect is the abrupt disappearance of the dominating large residue 
between 0.35 and 0.45 eV. At even larger excitations the cluster decays only into 
small fragments with masses less than 10. Looking at the thermodynamic tempera- 
ture T, (the right scale in fig. P, we see an S-shape in T(e*) around e* ~ 0.4 eV or 
T « 1200 K. This S-shape is connected with the disappearance of the large residue 
and is a typical signature of a microcanonical first order phase transition. 

The average masses Ml, M2 and M3 give of course only a rough indication 
about the real mass distribution at any given energy. Figure [| shows the detailed 
mass distributions for the same system as in fig. [1] for three individual excitation 
energies corresponding to just before, during, and after the phase transition. The 
picture does not show the monomers and does not distinguish between the fragments 
with the same mass but different charge. Before the phase transition (top panel) 
Na^j~ prefers to evaporate singly charged trimers and Na^~. This is due to the shell 
effects which we include for small clusters (N < 21). For lower excitation energies 
0.2 < e* < 0.25eV the Na^ channel is much more important than the Nag channel 
0. The predominance of the trimer channel is consistent with molecular-dynamics 



predictions for clusters with 4 < iV < 12, |TJ| and with the experimental findings 
in UTS] . In our case, however, the occurrence of Naf cannot be explained just by 
the least fission barrier, assumed in [15] . We do not restrict the number and type 



of the outcoming fragments; therefore, the choice of the system to fission instead 
evaporation or multifragmentation in the energy interval 0.2 < e* < 0.3 is not just 
an energetic, but mainly an entropic effect. 

The large residue in figure 0, top, has preferred masses Naj^, Na^ 7 and Na^, 
corresponding to the emission of exactly one fragment with mass N = 9, 3 or a 



neutral monomer. Inside the phase transition region (middle panel) we see a very 
broad distribution of residues combined with two dominating fragments, neutral 
Na4 and Na^~. The two available charges go mostly to the large residue and Na^~; 
other singly charged small fragments also occur but with much lower probabilities. 
One also finds channels with doubly charged residues with masses 47 to 61, but this 
decay channel is even less probable. 

After the phase transition (lower panel) no large residue is left and the system de- 
cays only into small fragments. The charges prefer to go to two Na<j~ clusters. Again, 
other singly charged clusters with masses from 3 to 15 exist, but their probabilities 
are at least one order of magnitude lower. 

Before studying Coulomb explosion we would like to test the importance of po- 
larization effects at different charges. Therefore we compare the conducting spheres 
(cs) calculation to the point charge (pc) calculation for Nayo" in figure [| A related 
question is that of the binding distance db between the clusters, because the polar- 
ization is strongly dependent on distances. In Appendix C we estimate dj, as 0.5 to 
1 A. In figure |3] we try both db values. We see that for Z = 2 the thermodynamic 
temperature, T(e*), is insensitive to db and polarization. At the same time one can 
see that the statistical uncertainty in T is about ± 70 K for e* < 0.4 eV. For larger 
energies the fluctuations become smaller. 

Let us now proceed to higher charges. Figure |] is the same as figure p] but 
for Na^. Note the change in temperature and energy scales. The disappearance 
of the large residue at e* « 0.09 eV coincides with a systematic small S-shape in 
T(e*), indicating a phase transition. For this charge, the signature of the phase 
transition is of the same order of magnitude as the statistical fluctuations of T(e*) 
in our calculations. For even higher charges the phase transition disappears. For 
Z = 8 the phase transition is shifted significantly towards lower excitation energies 
and temperatures for the simple reason that the total energy which is available 
to the fragmenting system is increased by the Coulomb self-energy Therefore the 
cluster needs much less additional excitation to undergo the phase transition. The 
lower transition temperature T indicates that by increasing the excitation energy, 



the additionally gained phase space is much less here than for Na^- This is quite 
natural, since the phase space for Z = 8 is more strongly restricted by the need to 
distribute all 8 charges among the small fragments. 

Looking in figure [| at the mass distributions for Na^J" before, during and after the 
phase transition, we see a slightly different picture from that for Na^. We find that 
at these three energies practically all the fragments, including the residue, are singly 
charged. In the top and middle panels, the small fragments are predominantly Najj~ 
and Na^, whereas, just after the phase transition (lower panel), they are Na^ and 
Nag". The charged trimers disappear because the system tries to put all available 
mass into charged fragments. At even higher energies the charged trimers appear 
again together with neutral fragments. 

To illustrate the role of polarization effects, we compare figure |5| {Z = 8) with 
figure [3] [Z = 2). The statistical fluctuations make the interpretation somewhat 
difficult, but two important things are still evident: First, looking at the solid and 
dotted curves for the conducting spheres interactions (cs) we see that the phase 
transition, present for df, — 1 A, is absent for df, = 0.5 A. It is intuitively clear, 
that there should be some charge Z which makes the cluster unstable even at e* = 
and it explodes into small fragments without any additional excitation. Taking the 
binding distance df, = 1 A this critical charge is Z = 9, as we show later. Re- 
ducing df, to 0.5 A we allow the fragments to be placed closer, so that attractive 
many-body polarization forces become even more important and cannot be com- 
pensated by repulsive forces any more. This additionally gained energy can be now 
spent to produce small fragments. The amount of this energy is in this case already 
sufficient to eliminate the phase transition. (Compare with figure [?[) Perform- 
ing the same calculation for db = 0.5 A for the point charge interaction (pc) and 
thereby neglecting the polarization effects, we still see evidence for a phase tran- 
sition (supported also by appearance and later disappearance of the large residue, 
as shown in figure ||). The second observation is that, outside the transition region 
(e* < 0.04 eV and e* > 0.09 eV), the temperature is systematically higher for con- 
ducting spheres than for point charges. For higher excitations, this shift is larger 



because of the appearance of neutral fragments, which are especially sensitive to at- 
traction by polarization. Taking polarization effects into account reduces the total 
Coulomb energy. The system then prefers to put this gained energy into the kinetic 
energy of the fragments, (See figure |9|, lower part.), increasing the temperature of 
the fragmenting configuration. 

Figure [5] shows the partitioning of the energy into Coulomb, kinetic, excitation 
and binding energy for decaying Na^J" and Na^. The energies are normalized with 
—Efot- We see that the Coulomb energy of the decaying configuration is negligible 
for Z — 2, but becomes absolutely dominant for Z — 8. Nevertheless, it is not the 
Coulomb effects but rather the change in mass partition that is significant for the 
fragmentation phase transition, and this is of the same type for both high and low 
charge. 

Extending now our results to a whole range of charges, we show the caloric 



equation of state T(e*) for charges Z=0...10 in figure 10. To help in orientating 
ourselves in this and the following diagram, we indicate some positions with A, B 
and C. The dotted curve (A,C) marks the fragmentation phase transition. High 
charges on the cluster suppress the transition signature and make it disappear for 
Z >8. 



To stress again the significance of the mass partition, figure |TT| shows the av- 
erage mass of the biggest fragment as a function of e* for charges Z=0...10. The 
dotted curve (A,C) again indicates the fragmentation phase transition. Since, in- 
dependently of the charge, the transition is signaled by the disappearance of the 
large residue we can identify the two phases simply by the presence or absence of 
one dominating large residue in the fragmented configuration. Since this is a micro- 
canonical first order phase transition, there is also a coexistence region in a certain 
energy and temperature range (S-shape in T(e*)). In this region, the fragmentation 
events can belong to one or both phases but outside it, only to one of them. 

The fundamental thermodynamic quantity is not the temperature, but (3. (See 
equations |].) Figure [L2| shows /3(e to t) in dependence on the total available energy, 
for Naf + for charges Z=2...10. We see that keeping the total energy e to t constant 
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and increasing the charge reduces the available kinetic energy and increases (3. This 
is a simple way to see how a non- energetic change of the system modifies the tem- 
perature. 

In this paper we studied the N-body phase space of Naf + -clusters as a function 
of charge Z and excitation energy. We conclude that in the statistical fragmentation 
process of metal clusters we can see micro canonical fragmentation phase transitions 
for different cluster charges. We can identify the two phases by the presence or 
absence of a big residue in the fragment configuration. Very high charges destroy 
the phase transition and force the cluster to explode into small fragments even at an 
excitation energy e* = 0. For Na 70 + this happens for Z > 8. The predicted charge 
and mass distributions in such a decay process might be also accessible experimen- 
tally. A question still open is whether a dynamical fragmentation follows the details 
of the phase space. This is connected to the question as to whether the dynamics 
is sufficiently ergodic or, in other words, whether an experiment is an equilibrated 
process. Different methods of exciting or charging the cluster might lead to differ- 
ent results concerning the equilibration condition. If this condition is reached, one 
should be able to see the signature of a phase transition in a suitable experiment. 

Appendix A: Polarization effects in many-body sys- 
tems and the approximation by the two-body inter- 
actions 

Coulomb interaction is a long-range interaction which is responsible for the overall 
correlation of the fragmenting system. This interaction contributes significantly to 
the total (conserved) energy. In our calculation we consider all the fragments to 
be conducting and spherical because alkali clusters are metallic even for very small 
sizes |Tj|. In a fragmenting many-body configuration the polarization effects are 



important. It is known how to treat the polarization of two spherical droplets from 
ITJI . Here we address the questions of how to treat polarization effects in a N-body 



system and how good the spherical approximation is for very small clusters. 

Generally, calculating Coulomb energy of a conducting sphere with charge Z and 
radius R in presence of a point charge Zp at distance d from its center we need to 
evaluate the interaction of an inhomogenious charge distribution on the surface of 
the sphere with this point charge Z P outside. Evaluating E cou i by the method of 
image charges we substitute this complicated interaction by the sum of two point 
charge interactions: first, Zp with the charge Zq in the center of the sphere, and 
second, Zp with its image charge Zi which is at distance dj from the center of the 
sphere fI7fl . The Coulomb energy resulting from the surface charge distribution is 
exactly equal to the E cou \ from the image charge method if 

Zi = —j z p, d i = -j-\ (3) 

Z C = Z-Z T , (4) 

because then the surface of the sphere is an equipotential surface. This is exactly 
the condition for the distribution of surface charge on a conducting sphere. 

In our case we want to calculate a Coulomb energy of several charged or neutral 
spherical clusters at close distances. We are going to iterate the image charge method 
which gives us at every iteration step a series of image charges which keep the surface 
of each conducting sphere at constant potential. The total interaction energy is the 
interaction of image charges of each sphere with all image charges in other spheres. 
Therefore the exact value of E cou \ can be achieved numerically with any desired 
accuracy by performing this iteration. 

Let us consider for simplicity first the case of only two clusters A and B with 
positive charges Za and Zb and radii Ra and Rp at a distance dAB of their centers, 



see figure [13|. Neglecting any polarization effects in the iteration step gives the 



Coulomb energy E as: 



^o = —, • (5) 

"AS 



In the first iteration step (1) we consider that the point charge Za induces a negative 
image charge Z Bl on the cluster B, and vice versa. The image charges are at 
distances o^i and dpi from the centers of clusters A and B and are on the line 
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connecting both centers. 

7 (i) _ Ra 7 , R\ . 7 (i) _ R b 7 R\ , , 

Z A1 — —-, Z B, «A1 — j , ^£1 — —j ^A, «B1 — j • yO) 

UAB UAB UAB UAB 

To account for charge conservation one needs to subtract the charge Z Bl in the 
center of cluster B and Z A1 in the center of cluster A. We get for the central charge: 

7(1) _ 7 7(1). 7(1) _ 7 7(1) ( 7 \ 

The Coulomb energy E\ after the first iteration step is then: 

7 (i) 7 (i) 7(1)7(1) 7(1)7(1) 7(1)7(1) 

dAB ~ dsi d-AB — d.A\ &AB ~ «A1 ~ «B1 d>AB 

In the second iteration step we need first to adjust the image charges at positions 
d>Ai and dsi, since the charges in the centers were changed. 

7 (2) _ Ra 7 (i). 7 (2) _ ^B 7 (i) rq v 

^Al — — ~> ^BOi ^Bl — —-, ^AO- V*) 

a AB dAB 

(2) (2) 

Further, the image charges Z A { and Z B { induce the second order image charges 

Z A2 an< l z m . 

7 ( 2 ) _ Ra 7 (2) , #j nn x 

Z ^2 - -j 1 Z B1> «A2 - "1 j , {i-V) 

a AB ~ a B\ a AB — «B1 



,-j, #g 7 (2) , _ ^B , 1 p 

Z B2 - _ 3 j — Z A1> «B2 - : — • (11- 



,(2) _ rx B (2 ) _ Jt B 
^1 ~A All «B2 — -. , 

«AB _ a A\ UAB — «A1 



For the Coulomb energy we get using cIao = g?bo — 0, 

2 7(2) 7(2) 

e 2 =y. , y* 3 ' • (12) 

ij=0 a AB — U-Ai — U-Bj 

In every next iteration step n we consider two aspects: 1) The charges in the 
center have been changed, therefore we need to adjust the image charges Z A { and 
Z B i, and 2) these image charges induce themselves second-order image charges Z A2 
and Z B2 , and so on up to Z A l and Z B ^. For the Coulomb energy we need to account 
for contributions from all the image charges in sphere A with all image charges in 
sphere B: 

n 7M y( n ) 

E n =z . M : Bi , • (13) 

j j =0 UAB — U-Ai — U-Bj 
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We can show (T^j that E n build a convergent sequence of energies. If charges Za 



and Zb have the same sign then 

E >E 2 >E 4 > ... E n ... >E 3 >E 1 . (14) 

Please note that the even subscripts are on the left side the odd ones on the right. 
The extension of the method to 3 or any number of conducting spheres is only a 
bookkeeping problem. Figure [P| shows, for example, the positions of image charges 



in the second iterations step for the case of three clusters. Let us concentrate 
on the left picture showing the complete collection of image charges. Besides the 
image charges on lines connecting the centers, which are represented by circles, we 
get image charges resulting from the interaction of the first-order images in one 
sphere with the other two spheres. The positions of these images are shown by 
triangles. Calculating energy we use equation [13| extending it to the general case by 
the summation over all the spheres. 

It is also clear that this method is valid if some of the conducting spheres are 
neutral. They are attracted due to polarization by the charged clusters. 

Our main aim is to find a representative set of configurations in phase space fl. 
For 1000000 configurations with 3 to 4 fragments it is no big problem to calculate 
the Coulomb interaction exactly. But for most of the excitation energies we need 
configurations with many more fragments. The exact treatment of the Coulomb 
energy makes the statistical task unsolvable within a reasonable computation time. 
Therefore we investigated whether we can separate the many-body interaction into 
the sum of two-body interactions. 



We first concentrate on a three-body case, see again figure [14]. Considering only 
the two-body interactions, right picture, we neglect the image charges symbolized 
by triangles in the left picture. 

We describe the positioning of three clusters by two variables, the distance dc 
of the center of the sphere C to the center of mass of A and B and cIab the distance 
between the centers of A and B. Figure [l^ shows the Coulomb energy surface for 
three Na^o clusters. The plateau at the small cIab and dc is a forbidden area where 



12 



the fragments would overlap. Figure 16 is the difference AE between the exact 



three-body interaction and the sum of the two-body interactions. The difference is 
smaller than 2%. Testing different combinations of masses N = 10 and N = 100 of 
the three singly or multiply charged fragments always yields AE < 5%. The same 
is true for one or two small neutral clusters in the presence of charged ones. This is 
well within the desired accuracy. 

Generally, the polarization part of the total Coulomb energy increases for larger 
clusters and smaller charges. The polarization part becomes especially dominant for 
a big neutral cluster in the presence of small charged ones at very short distances. In 
this case the approximation by a sum of two-body interactions gives an arbitrarily 
large deviation from the exact many-body calculation. In our fragmentation calcula- 
tion of a charged metallic cluster we never get configurations where the whole charge 
goes into small fragments aside one big neutral residue if the fragments have to be 
placed at surface distances larger than a binding distance db ~ 0.5A. Therefore we 
use the two-body conducting sphere approximation, evaluating the Coulomb energy 
in the MMMC calculation of metallic cluster fragmentation. 

Finally, to illustrate the role of polarization effects we have shown a calculation 
for the fragmentation of Na^j" and Na^ in figures ^| and |5|. We get systematically 
lower temperatures for the point charge interaction in comparison to the pairwise 
conducting spheres interaction. This means that taking the point charge interaction 
makes a smaller phase space accessible than in the case of conducting spheres. 

Appendix B: Polarization effects for dimers and 
trimers 

Most of the fragmentation configurations have many small fragments such as dimers 
and trimers. It is not obvious that treating dimers and trimers as spheres for cal- 
culating polarization effects is a good approximation. In order to clarify this, we 
calculated the interaction of the Na^ trimer-ion with a charged conducting sphere 
by means of the diatomics-in-molecules (DIM) method, modified in such a way as 
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to take into account the electrostatic interaction between the ion and the sphere as 



well as the mutual polarizabilities|2(|. The DIM method yields the potential energy 
surface of the trimer as a function of the inter-nuclear distances and can correctly 
describe the dissociation of the trimer-ion into Na^~ + Na or Na2 + Na + . In the 
presence of the sphere, the charge distribution in the molecule and the geometrical 
configuration may alter because of the electrostatic interactions. We found that the 
energy as a function of distance from the sphere is nearly the same as that calculated 
from the electrostatic energy alone provided that the charge distribution and geo- 
metrical configuration are allowed to relax freely. The energy also depends on the 
orientation of the trimer, there being two local minima at a fixed distance from the 
sphere: orientation A, an isosceles triangle having its apex between the sphere and 
the base, and orientation B, where the apex lies on the far side of the base as seen 
from the sphere. The energies for these two local minima bracket the electrostatic 
energy as calculated by treating the trimer as a sphere. Finally, when the trimer 
approaches too closely to the sphere, it becomes unstable because all of the charge 
flows onto the nearest atom, which is no longer bound to the neutral dimer. 



Figure [17] shows the two orientations A and B of a trimer corresponding to 
minima in the interaction energies. The curves show the dependence of the Coulomb 
interaction energy on the distance of the center of mass of the trimer to the surface 
of the cluster of radius 27 bohr and charge 5. The minimal orientation A is in nice 
overall agreement with molecular dynamics calculation |I3}| . The dot-dashed curve 
corresponds to the interaction between two conducting spheres, where the radius of 
the trimer-sphere was again taken as r tr imer = r s 3 1//3 . To put this approximation 



in relation to more simplified solutions, we show figure [TBI which is a blow-up of 



figure [L7|. The dotted curve is the approximation by point charge interaction for the 
charge localized in the centers, and the short dashed curve is when considering the 
large cluster as a polarizable sphere and the trimer as a point charge outside. 

This shows the validity of a spherical approximation for calculating the polar- 
ization energies for trimers. 

For a dimer the situation is even better, since there is only one minimum in 
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energy for the orientation radial to the cluster. In this case all the polarized charge 
moves as in the image charge method above. 

Of course, considering many-body configurations, the dimers and trimers cannot 
have the same orientation to all the fragments. Even though we cannot account for 
the individual orientations because of the demands on computer time, we know at 
least that the spherical approximation, which is like averaging over all orientations, 
gives a reasonable result. 

Appendix C: The minimal surface distances dj, be- 
tween the clusters 

Another much more important point is the minimal surface distance <4 for the 
clusters not to be bound to each other. It would be natural to assume that this 



distance is of the order of magnitude of the sum of the electronic spill-outs f21| , |22 
of the two interacting clusters. The problem is that spill-out is in general charge 
dependent and is probably strongly deformed by the presence of another charge 
outside the considered cluster. Therefore we obtained a rough estimate from the 
DIM calculation (see Figure |T7p, which also gives the forces experienced by the 
three trimer atoms. The calculation break down when the attraction of one of the 
trimer atoms to the cluster is larger than its binding within the trimer. This minimal 
distance gives an estimate for the lower limit of d b as d b « 0.5 to lA. 
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Figure 1: Na^: Average three largest masses Ml, M2, M3 (left scale) and the 
thermodynamic temperature (right scale) as a function of specific excitation energy 
e*. df, = I A, pairwise conducting spheres interaction. 

Figure 2: Na^j": Mass distributions for different specific excitation energies e* 

Figure 3: Na^: Thermodynamic temperature T(e*) for different db values and 
conducting spheres (cs) versus point charge (pc) Coulomb interactions. 

Figure 4: Like figure [I] for Na^. 

Figure 5: Like figure |2|for NayJ". 

Figure 6: Like figure ||| for Na^J"- 

Figure 7: Like figure ^ for d\, = 0.5A. 

Figure 8: Like figure |^for point charge interactions. 

Figure 9: Partition of the total available energy for fragmenting Na^j~ and Na^j". 

Figure 10: Na 70 + : Caloric equation of state T(e*) for charges Z=0...10. The dotted 
curve indicates the fragmentation phase transition. 

Figure 11: Naf + : Average mass of the biggest fragment in dependence on e* for 
charges Z=0...10. The dotted curve indicates the fragmentation phase transition. 
The Points A, B and C indicate the corresponding places in figure |l~0| . 

Figure 12: Naf + : (3{e*) for charges Z=2...10. 

Figure 13: First iteration steps for evaluation of Coulomb interaction energy by the 
method of image charges in case of two conducting spheres A and B with charges 
Za and Zg. 

18 



Figure 14: The second iteration step for evaluation of Coulomb interaction energy 
by the method of image charges for three conducting spheres A, B and C with 
charges Za, Z b and Z C - 



Figure 15: Exact three-body Coulomb energy E of three spherical conducting 
clusters with mass N=10 and charge Z=l as a function of different relative positions, 
see text, Appendix A. 



Figure 16: Difference AE between E, figure [15] and the sum of the two-body 
interactions of three spherical conducting clusters with mass N=10 and charge Z=l 
as a function of different relative positions, see text, Appendix A. 



Figure 17: Coulomb energy for a charged sodium trimer in a Coulomb field of a 
spherical cluster with radius 27 bohr, Z=5 for two orientations A and B, see text, 
Appendix B. The dot-dashed curve is a spherical approximation of a trimer. 



Figure 18: The three curves from figure |TT|are shown together with two other 
approximations for Coulomb energy: Dotted curve considers both, a sphere and a 
trimer having fixed point charges in their centers. Dashed curve shows the point 
charge approximation just for the trimer in presence of conducting sphere. 
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